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I. INTRODUCTION 



A. BACKGROUND 

I. Tolstoy [Ref. 1, Ref. 2], and A. Tolstoy and 1. Tolstoy [Ref. 3] have 
recently analyzed the multiple-scattering of sound from a collection of 
compact (ka « 1) monopole scatterers (e.g. gas-filled bubbles, gas- 
filled balloons, and thin shells) in water, insonified by a simple 
harmonic plane-wave at frequencies close to the resonance frequency 
of an individual scatterer. The response of these systems was 
quantified by computing the ratio of the scattered pressure on the 
surface of a scatterer in the presence of the other scatterers to that of 
the incident plane wave. This ratio was defined as the amplification 
factor, |x. They investigated the value of this amplification factor for 
various simple geometrical arrangements of monopole resonant 
scatterers, and showed that amplification factors of the order of 10^ 
were to be expected for air bubbles in water at one atmosphere 
ambient pressure (The amplification factor of a single bubble is about 
70.). These systems were said to be "quasiresonant". 

1. Tolstoy (Ref. 4] has also recently analyzed the strength and 
directionality of the resonant multiple-scattering of sound from a 
linear array of compact monopole resonant scatterers insonified by a 
simple harmonic plane wave. He reports that "quasiresonance" effects 
exhibit strong directionality in this system, with amplification factors 
of up to 700. The potential application of resonant multiple- 
scattering, or "quasiresonance", in Tolstoy’s words, to improve the 
performance of a low frequency active sonar transducer is the subject 
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of the research described in this report. 



B. OBJECTIVES 

The objectives of this research are multiple. The first is to 
reproduce I. Tolstoy's results using network analysis, a tool familiar to 
the transducer and array designer. It will be shown that network 
analysis gives results identical to Tolstoy's analysis. The existence of 
quasiresonance will then applied to a circular array of scatterers, 
where the incident plane wave is replaced by an active transducer 
located on the axis of the circle and radiating spherical waves. The 
acoustic advantage gained in this system will be quantified by the gain 
in the radiation resistance seen by the trcmsducer by the introduction 
of the array of scatterers. An increase in radiation resistance enables 
the transfer of more power to the acoustic field for a given surface 
velocity of the transducer (High power low frequency active sonar 
transducers tend to be resonant devices that are displacement 
limited.). The cases of the transducer in the plane of the circle and 
displaced by one-quarter wavelength out of the plane will be analyzed 
for both their increase in radiation resistance and gain ir directivity 
over that of a single transducer. 
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II. THEORY 



A. INTRODUCTION 

For simple configurations of compact resonant scatterers insonifled 
by a simple harmonic plane wave, Tolstoy predicted the existence of 
so-called "quasiresonance" using multiple-scattering theory 
(Quasiresonance is Tolstoy's term for the dramatic enhancement of 
the scattered field of a resonant scatterer as a result of resonant 
multiple-scattering.). Tolstoy's analysis will be reviewed and applied 
to the simplest case of the scattering of a plane-wave axially incident 
upon a pair of resonant scatterers in Section B. In Section C the same 
problem is treated using network analysis. It will be shown that this 
method produces results identical to Tolstoy's. 

B. TOLSTOY'S ANALYSIS OF RESONANT MULTIPLE- SCATTERING 

Tolstoy's analysis of the multiple-scattering of a plane wave incident 
upon a collection of compact resonant scatterers proceeds using the 
linear theory of acoustics. The following development parallels that of 
Reference [1] and Reference [3], with the difference being that here 
all time-harmonic functions are taken to be of the form e)“k This will 
result in minor differences in the appearance of otherwise equivalent 
equations. 

Consider first the case of a single resonant scatterer insonified by 
an incident plane-wave. The scatterer is taken to be an air-filled 
bubble of radius a, which is assumed small compared to the acoustic 
wavelength (ka « 1). From Klnsler, Frey, Coppens, and Sanders [Ref. 
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5. p. 230], the ratio of the scattered pressure on the surface of the 
bubble to the free-field pressure of an incident plane-wave is. in the 
notation adopted here (note: the minus sign in Eq. 2-1, which is not 

present in Eq. 10-22 of Ref. 5, results from the assumed convention of 
positive velocity directed into the water); 



where: 

pifngle(a) = scattered acoustic pressure on the bubble surface. 



pf = incident free field acoustic pressure, 

Rr = bubble mechanical radiation resistance, 
mr = bubble mechanical radiation mass, 

Rm = bubble internal mechcuiical resistance, 
s = bubble internal mechanical stiffness. 



The scattered pressure at any distance r, plfngle^*")* related to the 
scattered pressure at the surface by 



where the approximation ka « 1 has been used. Combining Eq. 2-1 
and Eq. 2-2 results in the following expression for the ratio of the 
scattered acoustic pressure at some radius r from the scatterer, 
Pl?ngle(r). to that of the incident pressure, p^: 



Pifngle(^) 



(Rm + F^) + j(comr - s/co) ’ 



Rr + Jcomr 



( 2 - 1 ) 



Pffngle^^^ _ a ~ 

Pflngle(a) ^ e-Jkr “ jkr (1 -Jka) ’ 



( 2 - 2 ) 




(2-3) 



where: 



. ^ (c0§/(0^ -l) 
ka 
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and 



top ^ 
to 



3yPq 

pc2 (ka)^ 



(2-4) 



In deriving Eq. 2-3, the following assumptions and substitutions have 
been made: 

(1) Surface tension is neglected, so the only mechanical 
stiffness of the bubble is due to internal gas (air) pressure: 

(2) The bubble internal mechanical losses due to viscosity and 
thermal conductivity are neglected compared to the acoustic 
radiation loss, i.e. Rm « Rr: 

(3) Rr = (Oomr/Q: 

(4) 9 = l/(ka); 

(5) S = too^mr: 

(6) ka « 1: 

(7) to = top. 



It follows from Eq. 2-3 that the ratio of the scattered acoustic 
pressure on the single scatterer's surface to that of the incident free- 
field pressure can be written: 



^ _J_ I 1 \ e-jka 

p' Jkall-jd 



(2-5) 



Tolstoy introduces the quantities Fa and A by rewriting Equation 2- 
3 as Equation 2-6. 



Fa 





( 2 - 6 ) 



where: 

A = -L. 

1-jC 



Fg is the normalized scattered acoustic pressure for a single scatterer 
insonified by an incident plane-wave and A is the single-scattering 
amplitude coefficient. 
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Now consider a system of N identical interacting scatterers 
insonifled by a plane-wave. Let pfip(r) represent the scattered acoustic 
pressure. An equation for p?f(r)/pf equivalent to that of Equation 2-6 
for a single scatterer, can be written as Equation 2-7 for the multiple- 
scattering case: 

= (2-7) 

pfi Jkr 

is the normalized scattered acoustic pressure from scatterer n, 
insonified by both the incident plane-wave and the fields of the other 
scatterers, p|i is the free-field pressure of the incident plane-wave at 
scatterer n, and is the multiple scattering amplitude coefficient for 
scatterer n. 

The response of a scatterer to multiple incident pressure fields can 
be evaluated by summing up the response to each individual incident 
pressure field [Ref. 3]. For a system of N scatterers, this results in 
Equation 2-8 for the scattered field at the surface of the n^L resonator: 

B„ = A„ / 1 -h X Bn. , (2-8) 

J^n jKinm j J^an 

where: Vnm the phase difference in the incident plane wave 

between resonator m and n, and m = 1...N. The subscript notation for 
A has been added to allow for a system of different- sized scatterers. 

Tolstoy defines the amplification factor |i (different from the 
scattering amplitude coefficient mentioned previously) to be the 
magnitude of the ratio of the scattered acoustic pressure on the 
surface of a scatterer, pftc^ ^he free-field acoustic pressure of the 
incident plane-wave, pfi. With this definition. Equation 2-9 follows 
from Equation 2-7: 
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M^n - 



(2-9) 



pfP\ _ |Bnl 

pfi kan 

Tolstoy denotes the amplification factor for a single scatterer as go- 
Equation 2-10, then, follows from Equation 2-6. 







_ l^nl 

kan 



( 2 - 10 ) 



From Equations 2-9 and 2-10, the gain In the amplification factor 
due to multiple-scattering, |i/|Xo. is expressed In terms of Bn/An by 
Equation 2-11: 



gn 

l^On 



P?T 

Pi&iglen 



= |Bn/Anl . 



( 2 - 11 ) 



B.l. APPLICATION OF TOLSTOY'S METHOD TO A SYSTEM OF 
TWO RESONANT-SCATTERERS 

As an example application, Tolstoy considered a system of two 
Identical resonators Insonlfied by a plane-wave traveling along the axis 
joining their centers, as shown in Figure 2-1. 

Resonator of 




Figure 2-1. A two-resonator system insonlfied by a plane-wave 
along the axis of the resonators. 

Application of Equation 2-8 to the situation shown in Figure 2-1 
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results in the set of Equations 2-12: 



B eilhi = A + A B 2 e-Jkl . (2-12) 

jka jka jka ^ jkl 

B2 = A + A Bi e-Jkl . 

jka jka jka jkl 

For a system of identiCcil scatterers, the values of Ai and A 2 are equal; 
therefore, the subscript on A in Equation 2-12 has been omitted. 
Since all variables in Equation 2-12 are known except Bi amd B 2 , one 
can solve for Bi/A in terms of A amd kl, resulting in Equation 2-13: 



Bi/A = 



1 + (a e-Jkl 



(2-13) 



Substituting the expressions for ^ and cd/ooq for an eiir-filled bubble into 
the expression for A, Bi/A is given in terms of the variables ka, kl, c, 
Pq, and p by Equation 2-14: 



1 - e-jW 



Bi/A = 



_ Pf _ 





ka ej(ka - kl) 
kl 




^pc 1 


|(ka)^ - 1 - j 


fka- 1 


1 


^ pc2(ka) I 



Pffngle 



1 - 





ka e)(ka - kl) 
kl 


\2 


1 


|(ka)^ + j| 


fka- 1 


I) 


i pc2(ka) 1 



(2-14) 



For bubbles just under the surface of the ocean: 

Po = 10132 Pa (1 atmosphere) 
p = 998 kg m~3 
c = 1498 m s~i 

As observed by Tolstoy [Ref. 3), plots of Equation 2-14 display 
spurious peaks. Following Tolstoy, then, the data used to produce the 
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graphs of IBi/AI versus ka and kl, plotted in Figures 2-2a,b, have been 
treated to limit peaks in bounded intervals of width less than 2ka 
(According to Tolstoy, peaks narrower than 2ka, where 2a equals the 
resonator diameter, are not meaningful, since the pressure over a 
compact resonator is approximately constant.). The process used 
involved comparing three consecutive values of IBi/AI in the kl 
direction for each value of ka (Three consecutive points in the kl 
direction span approximately 2ka.). If the value of IBi/AI for the 
middle point was greater than either of its neighboring points (a peak) 
then this value was set equal to the mean of the two neighboring 
points. This procedure was repeated several times to obtain the 
desired smoothing. 

Figures 2-2a,b show a maximum response of IBi/AI equal to 5.5 at 
approximately ka = 0.01393 and kl = 0.60. As a result of multiple- 
scattering, then, the scattered pressure amplitude of the first bubble 
is 5.5 times greater in the presence of the second bubble than in its 
absence. Additionally, the value of ka at the maximum response is 
shifted from that for a single resonant bubble of 0.01379. 
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Figure 2-2a,b. Surface (a), and contour (b) plots showing iBi/Al 
for the double-resonator system as a function of ka and kl for 
axial Incidence of the plane-wave. 
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C. NETWORK ANALYSIS OF RESONANT MULTIPLE-SCATTERING 

Figure 2-3 displays a network model of N interacting scatterers 
insonified by an incident acoustic wave (Note: pressure (p) and 

volume velocity (Un) are chosen as the mechanical variables, and the 
direction of positive velocity is clockwise.). 



Resonators Acoustic Field 




Figure 2-3. Network model of N interacting resonators 
insonified by a plane-wave. 



The network equations for Figure 2-3 are: 

N 

Pn = pfi + ^ ZnrnUm Pn = “ Un , 



m = 1 



(2-15) 



where: 

Pn = Surface pressure at element n, 
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pf, = Free field pressure at element n due to an 
externally incident field only, 

Un = Volume velocity of element n, 

Zgn = Acoustic impedance of element n, 

^nn = Acoustic radiation impedance seen by element 

n with all other elements blocked (Um ^ n = 0: m=l...N), 
Znm = Acoustic transfer Impedance from element m to 
element n. 

Neglecting the effects of surface tension, thermal relaxation, etc., 
the intrinsic acoustic impedance of a small (ka « 1) gas-filled bubble 
can be approximated by one of reactance only, given by Kinsler, Frey, 
Coppens and Sanders [Ref. 5) as Equation 2-16: 



where: 

a„ = radius of element n; 
c = speed of sound in water; 
k = wave number, co/c; 

Pq = ambient pressure of environment; 

Sn = surface area of element n; 

y = ratio of specific heats (air). 

The self acoustic radiation impedance seen by element n, Znn. can 
be approximated by Equation 2-17. 

= pc ((kanP+jkan) S“^ , (2-17) 

where p is the density of water. Equation 2-17 is the acoustic 
radiation impedance presented to a compact source with no other 
active or scattering elements present. It is also the acoustic radiation 
impedance seen by a compact source in the presence of other 
compact sources when these sources are blocked (Ujn n = 0: = 
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1...N). 



The acoustic transfer impedance from element m to element n is 
given by Equation 2-18: 



- ^nm) . (2-18) 

klnm 



In matrix form, the network equations (Equation 2-15) can be 
written as Equation 2-19: 

{Pn* = [Znm] (Um) + |pfil ’ ^^i.n = 1...N , (2- 19) 



where: 

IXjil: n = 1...N A IxN column matrix, N elements, 

[Xnm]« n,m = 1...N A NxN square matrix, elements. 

Substituting Pn = ~ Z^Un and rearranging results in Equation 2-20: 

1“ Pn! = [Znm + Zn5nm] (Um) I m,n = 1...N , (2-20) 

where is the Kronecker delta function: 



mm 



J 



1; m = n I 
0; m n j 



C.l. APPLICATION OF NETWORK ANALYSIS TO A SYSTEM OF 
TWO RESONANT-SCATTERERS 

Modeling the system of two scatterers shown in Figure 2-1 as a 
network. Equation 2-20 becomes Equation 2-21: 



-P^ 




zf5 + zf 


zf5 


Ui 


- p£ 




z|J 




U2 



Using the relation pftc = ZnnUn to eliminate Un. the above equations 
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can be written in terms of the scattered acoustic pressure at each 
individual resonator, pgp, and the incident plane-wave pressures, as 
given by Equation set 2-22: 



pf = (p^ + Ci2 pr ) Dl^ . 
pr = (p£ + C21 pf )D2^ , 



( 2 - 22 ) 



where 




<73 r 
■^onm 



and 




For a system of identical resonators, Cnm and Dn are the same for all m 
and n; therefore the subscript notation is not required. Solving 
Equation 2-22 for the scattered acoustic pressure at each resonator in 
terms of the incident plane-wave pressure results in Equation 2-23 for 
the scattered pressure at resonator one: 



The scattered acoustic pressure at the first resonator in the 
absence of the second resonator, pffnglei’ obtained using Equation 
2-23 by letting the transfer coefficient C equal zero, resulting in 
Equation 2-24: 



Dividing equation Equation 2-23 by Equation 2-24 results in 
Equation 2-25 for the ratio B\/A as defined by Tolstoy: 




(2-23) 



Pslnglei - ^ 



P{ 



(2-24) 



D 
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(2-25) 



Bi/A = 



- Pf - 



1 + (p^M) 



Pffngle] 



-(gr 



Inserting Equations 2-16 thru 2-18 into Equation 2-25, one obtains 
an equation for Bi/A in terms of the variables ka, kl, c, Pq, p^, p^, and 

P, as given by Equation 2-26: 



1 "(pS/pi) 



Bi/A = 



- pf _ 



ka^JIka - kl) 
Id 






SC 

Psinglei 



ka^(ka - kl) 
Id 



. (2-26) 



pc((kaP+j(ka--^jj 



The incident plane wave pressures p^ and p^ are of equal magnitude 
and differ only by the phase factor d*^i 2 , resulting in Equation 2-27: 



= e-jkl . (2-27) 

P\ 

Substituting this expression into Equation 2-26, one obtains Equation 
2-28: 



1 - e--Jkl 



b./a = ^^ 

Psinglei 



ka ^(ka - kl) 
Id 



pc 






1 - 



ka ej(ka - kl) 
kl 



(2-28) 



pc 



((kaP 



/ 



+ j 



ka - 



3yPo 

pc2(ka) 



Equation 2-28 is identical to Equation 2-14, demonstrating the 
equivalence of Tolstoy's method and network analysis. 
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D. NETWORK ANALYSIS OF A TRANSDUCER INTERACTING WITH 
MULTIPLE RESONANT-SCATTERERS 

Network analysis will now be applied to the problem of a compact 
transducer in the presence of an array of interacting, compact 
resonant scatterers. The procedure is the same here as in section C, 
except that the incident plane wave is replaced by an outward- 
traveling spherical wave whose source is a compact transducer. A 
network model of the transducer and N resonator system is shown in 
Figure 2-4. 



Acoustic Field 



Transducer 



Resonators 




Figure 2-4. Network model of a single transducer and N 
interacting resonators. 

The netv'ork equations for Figure 2-4 are: 
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po=TI-Z^°Uo. Pn = -Z^Un; n = 1..N, (2-29a) 

Pn “ ^ I m,n — 0..N, (2-29b) 

m 



where 

T = transduction coefficient: 

= transducer open circuit (I = 0) acoustic impedance, 



and the subscript n = 0 denotes the transducer, n 0 a resonator. 

The transducer open-circuit acoustic impedance (Z®°) and 
transduction coefficient (T) are shown only for completeness and will 
not be considered further, as the quantity of interest will be the 
radiation Impedance presented to the transducer. 

Using the network Equations 2-29, the matrix equation for the 
network model shown in Figure 2-4 is given by Equation 2-30: 



po 




Z85 IzgTor 


(Pnl 




(zS&l [zf,y 



Uo 

(Uml 



n, m = 1...N 

> 



(2-30) 



where 



IXj^l, n = 1...N denotes a IxN column matrix, N elements, 
[Xnml. = 1...N denotes a NxN square matrix, elements. 



and the superscript T denotes transpose. Replacing the 1 x N surface 
pressure matrix, {pnl. with (- ZnUnI, and then normalizing by Uq, the 
transducer volume velocity. Equation 2-30 becomes Equation 2-31: 



Zfo^tal 




(0) 





m 

|zgb| 



where 



„1 



(Tireh 
|Um I 



n.m = l...N.(2-31) 



Zfotal = total acoustic radiation impedance seen by the 
transducer; 
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Um^ = ratio of element m’s volume velocity to that of the 
transducer (Um^ = Um/Uo). 



Equation 2-31 can be written as two separate matrix equations, 
Equations 2-32a,b: 

Zf^tal = Z66 + . (2-32a) 

(0) = {^5j + [Znm + Zn5nm] {Um^} • (2-32b) 

Solving Equation 2-32b for {Um^l by matrix inversion and substituting 
the result into Equation 2-32a, one obtains Equation 2-33, the total 
acoustic radiation impedance seen by the transducer. 

Zfotal “ Z§5 ~ {Zmo}^ [Znm + Z^^ntn] jZnb} • (2-33) 

Equation 2-33 expresses this impedance in terms of the radiation and 
transfer impedances of the individual elements. 

The total acoustic power radiated by the transducer is given by 
Equation 2-34. 



n = i|Ud^Re(Zfo"tal). (2-34) 

The gain in power delivered to the acoustic field due to the addition of 
the sympathetic resonators equals the ratio of total power delivered 
with the resonators present to the power delivered when they are 
absent, Ilo = ;^ | UqI^ Re(z56 )• For a displacement-limited source, 
typical of a sonar projector operated near its resonant frequency, this 
ratio reduces to the ratio of the radiation resistances for the two cases. 
This gain factor is expressed as Equation 2-35: 



Gd = F.e(Zfotal ) 

Re(Z65) ■ 



(2-35) 
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III. ANALYSIS 



A. INTRODUCTION 

The theory developed in the previous chapter will now be applied 
to two specific geometric configurations of a transducer and N 
identical resonators: a "planar" and a "conical" configuration. The 

resonators will be treated as simple gas-filled bubbles. It is assumed 
throughout that the acoustic wavelength is much larger than the 
dimensions of either the transducer or the resonators. For initial 
calculations, ambient pressure is taken to be one atmosphere, and the 
radius of the transducer is taken to be equal to that of the resonators 
(kao = kan). Later, calculations will be presented for the case of an 
ambient pressure of 50 atmospheres, and the restriction kao = kan will 
be relaxed. 

B. ANALYSIS OF PLANAR CONFIGURATION 

The geometry of the first arrangement to be analyzed is shown in 
Figure 3- 1 . It consists of N identical resonators, equally spaced 
around a circle of radius InO- This geometric arrangement will be 
referred to as the planar configuration. 

For the configuration shown in Figure 3-1, Equation 2-35 becomes 
Equation 3-1: 

~ [Znm + Z^^nni] 4 (3-1) 

Re( Z85) 



27 





Figure 3-1. Resonators equally spaced around transducer 
in a planar configuration of radius InO. N=4. 



Using Equations 2-16, 17, and 18, and canceling the common surface 
areas Sn, Equation 3-1 can be written in terms of specific acoustic 
impedances as Equation 3-2: 



P ^ Re(zg5 - (zMo^ + zaSnm] ^ . 

Re( z86) 



n,m = 1...4 



(3-2) 



where: 



^86 = pc ((ka)^ +Jka) , 

+ zft5nm = pc ((ka)^ + jka) - j ; n = m, 

c(ka) 

zMn + zftSnm = zftlfn = pc ({kap + Jka) d(ka - klnm) ; n m, 

z^o = pc ((kap + jka) el Ika - klno) 

klnO 



and Inm is the distance between element n and m, the subscript 0 
denoting the transducer. Using simple trigonometry, klnm for n,m 
0 can be written in terms of m, n, and klno for the planar case as 
Equation 3-3: 
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(3-3) 



klnm = klnO 2 sinj^ jn - m|j ; n,m = 1...N. 

Using Equations 3-2 and 3-3, the gain in radiation resistance seen by 
the transducer, Gr, can be expressed in terms of ka and klno- 

The gain in radiation resistance for the N = 4 planar configuration 
can be plotted in the form of surface and contour plots using Equation 
3-1, as shown in Figures 3-2a and 2b. It should be pointed out that 
here and in all subsequent surface plots of Gr, no spurious peaks were 
observed, and no smoothing has been applied. Inspection of these 
graphs reveals two local maxima, one of magnitude 4.34 at ka = 0.975 
kar, klnO = 0-3, and another of magnitude 1.79 at ka = 1.005 kar, klno 
= 2.7, where kar equals 0.01379, the resonance value of ka for a single 
air-filled bubble at one atmosphere ambient pressure. The local 
maximum at ka = 0.975 kar, klnO = 0-5 is an artifact of the range of 
klnO chosen. In fact, the gain in radiation resistance diverges as 
(klno) ^ as klnO approaches zero. This case is uninteresting and will 
not be considered further. The interesting local maximum at ka = 
1.005 kar, klnO = 2.7 represents a gain of almost a factor of two in the 
radiation resistance seen by the transducer. Note that kl = 
corresponds to a separation distance of one-half wavelength. For the 
N = 4 planar configuration, then, the maximum gain in radiation 
resistance is obtained when the resonators are located a distance just 
under one-half wavelength from the transducer. 

The gain in radiation resistance for the planar configuration where 
N is set equal to eight and fifteen is shown in Figures 3-2c thru 3-2f. 
These plots display similar characteristics to the N = 4 case. All show 
a local maximum of magnitude approximately two at ka = kar = 
0.01379 and klnO = 
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The local maximum values for Gr at klnO s ti, ka = kar = 0.01379 for 
various values of N ranging from 2 to 25 and their associated 
parameters are tabulated in Table 3-1. Two major points can be 
observed from Table 3-1. First, the resonant bubbles' volume velocity 
leads the transducer's by a phase shift of 90’’ (equivalent to - 270°). 
This phase shift is created in two parts. Since the separation between 
transducer and resonator is approximately one-half wavelength, this 
introduces a shift of - 180° in phase due to simple wave propagation. 
The remaining - 90° is introduced at the surface of a compact 
resonating bubble between the incident and scattered wave, as 
predicted by Equation 2-24. Second, note that little is gained as one 
increases the number of resonators beyond about four. 



N 


Gr 




klfxiax 


mag(U{f‘) 


arg(U{f‘) 


2 


1.26 


0.01382 


2.65 


0.395 


90.0° 


4 


1.87 


0.01383 


2.82 


0.649 


90.0° 


6 


1.96 


0.01385 


2.87 


0.474 


90.0° 


8 


1.96 


0.01386 


2.87 


0.356 


90.0° 


10 


1.96 


0.01386 


2.87 


0.285 


90.0° 


15 


1.96 


0.01385 


2.87 


0.190 


90.0° 


20 


1.96 


0.01381 


2.87 


0.142 


90.0° 


25 


1.96 


0.01376 


2.87 


0.114 


90.0° 



Table 3-1. Values of Gr for N ranging from 2 to 25 and the 
parameters associated with the local maximum at klno = tc, 
ka = kar = 0.01379. Planar configuration. 
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a. 





Figure 3-2a,b. Surface and contour plots of Gr as a 
function of ka and klnO for N = 4. Planar configuration. Pq 
= 1 atm. 
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o 




c. 




Figure 3-2c.d. Surface and contour plots of Gr as a 
function of ka and klno for N = 8. Planar configuration. Pq 
= 1 atm. 
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Figure 3-2e,f. Surface and contour plots of Gr as a 
function of ka and kInO for N = 15. Planar configuration. 
Pq = 1 atm. 
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B.l. PLANAR CONFIGURATION BEAM PATTERNS 

Polar plots of the beam patterns can be developed for the planar 
configuration using the axis system defined in Figure 3-3. The origin 
is located at the transducer, with the resonators located in the x,y 
plane and with 0 being the angle measured from the positive y axis in 
the x,y plane. 

Figure 3-4 clearly demonstrates the acoustic advantage of the 
planar configuration over that of a lone transducer. The zero dB 
reference is taken as the field of an identical transducer with equal 
volume velocity, in the absence of the resonators. 



A z axis 




y axis 



Z 



Z 




X 



y.z plane 



z 0 =45 °plane 



Figure 3-3. Axis system for beam patterns plots. Planar 

configuration. Note the definitions of the y,z and z,0 = 45° 
planes. 
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z axis 



o 




Figure 3-4. Beam patterns in the y.z and z,0 = 45'* planes. 
Planar configuration. Pq = latm. Scale in decibels (dB). 



C. ANALYSIS OF CONICAL CONFIGURATION 

The second geometric arrangement analyzed was chosen to take 
advantage of the 90'* phase shift between the volume velocity of the 
resonators and the transducer for maximum radiated power. If the 
transducer is offset from the plane of the resonators by one-quarter 
wavelength along the positive z axis, the radiation from the resonators 
will constructively Interfere with that of the transducer along the 
positive z axis. This configuration will increase the directivity of the 
system. This arrangement will be referred to as the conical 
configuration, and is shown in Figure 3-5. 



35 




Resonator n 



Figure 3-5. Resonators equally spaced around transducer 
in a conical configuration at radius R. N = 4. The 
transducer is displaced out of the plane of the resonators 

by a distance of X/4. 

The value of Gr can again be expressed in terms of ka and kl^o 
using Equation 3-1, where Equation 3-4 is used to express klnm in 
terms of n, m, and klnO for fho conical configuration. 



The gain in radiation resistcince for the conical configuration for the 
cases N = 4, 8, and 15 are plotted in the form of surface and contour 
plots, as shown in Figures 3-6a through 6f. Inspection of these plots 
reveals the existence of the same local maximum as the planar 
configuration of magnitude slightly less than 2 at klno = and ka = 
0.01379. Note that, unlike the planar case, the kl axis of Figures 3-6a 
through f begins at kl^o = 1-75 vice 0.5. This is due to the 
displacement of the transducer in the positive z direction. This 
displacement requires kl^o be at least n/2. 

The local maximum values for Gr at klnO s n, ka = 0.01379 for 




n,m = 1...N. (3-4) 
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various values of N ranging from 2 to 25 and their associated 
parameters are tabulated in Table 3-2 for the conical configuration. In 
each case the value of Gr for the conical configuration is slightly less 
than for the planar configuration. Also, the value of klmax Is slightly 
greater for the the conical configuration, but is still approximately 
equal to n (A value of n for klno corresponds to a value of 30 degrees 
for the angle a in Figure 3-5.). 



N 


Gr 






mag(U|r‘) 


arg(U-‘) 


2 


1.29 


0.01382 


2.77 


0.429 


90.0° 


4 


1.78 


0.01383 


3.08 


0.601 


90.0° 


6 


1.82 


0.01385 


3.12 


0.428 


90.0° 


8 


1.82 


0.01386 


3.12 


0.321 


90.0° 


10 


1.82 


0.01385 


3.12 


0.257 


90.0° 


15 


1.82 


0.01383 


3.12 


0.171 


90.0° 


20 


1.82 


0.01378 


3.12 


0.128 


90.0° 


25 


1.82 


0.01372 


3.12 


0.103 


90.0° 



Table 3-2. Values of Gr for N ranging from 2 to 25 and the 
parameters associated with the local maximum at kl s tc, 
ka = klr = 0.01379. Conical configuration. 
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Figure 3-6a,b. Surface and contour plots of Gr as a 
function of ka and klno for N = 4. Conical configuration. 
Pq = 1 atm. 
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c. 




Figure 3-6c,d. Surface and contour plots of Gr as a 
function of ka and klnO for N = 8. Conical configuration. 
Pq = 1 atm. 
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n 





Figure 3-6e,f. Surface and contour plots of Gr as a 
function of ka and kl^o for N = 15. Conical configuration. 
Pq = 1 atm. 
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C.l. CONICAL CONFIGURATION BEAM PATTERNS 

To demonstrate the improved directivity of the conical 
configuration over that of the planar system, polar plots can be 
developed using the axis system as defined in Figure 3-7, where the 
transducer is displaced in the positive z direction from the plane of 
the resonators. 




Figure 3-7. Axis system for beam pattern plots. Conical 
configuration. 

Figure 3-8a shows the beam pattern of the conical configuration in 
both the y,z and the z,0 = 45° planes. As expected, the symmetry in 
the y,z plane about the z axis is lost, with the sound pressure level 
(SPL) in the positive z direction being 8 dB greater than in the 
negative z direction. If one compares the beam patterns of the conical 
configuration to that of the planar configuration, shown in Figure 3-8b, 
one observes that the SPL in the positive z direction for the conical 
configuration is 1 dB greater than for the planar configuration, and 10 
dB greater than for one transducer alone. Figure 3-8b clearly shows 
that significant directional discrimination can be achieved with only a 
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minor decrease in total power gain, Gr. 



z a^is 




Figure 3-8a. Beam patterns in the y,z and z,0 = 45" planes. 
Conical configuration. Pq = 1 atm. Scale in decibels (dB). 
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z axis 




Figure 3-8b. Beam-patterns in the y.z plane for the planar 
and conical configurations. Pq = 1 atm. Scale in decibels 
(dB). 



D. OPTIMUM CONFIGURATION 

The results of all calculations performed have been presented in 
terms of the dimensionless radius and spacing, ka and kl„o. The 
actual dimensions for the radius a and the separation distance Imax for 
a given frequency are obtained using Equation 3-5: 



T _ kac 



InO = 



hlnpC 

27tf 



(3-5) 



Taking ka = kar = 0.01379 and klnO s in Equation 3-5 for the 
optimum values of a typical system at an ambient pressure of one 
atmosphere, a plot of the resulting dimensions a and I^q ss a function 
of frequency can be made, as shown in Figures 3-9a and 3-9b: 
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10 




Figure 3-9a. Optimum resonator diameter (a) at resonance 
versus frequency. Pq = 1 atm. 




frequency (hz) 



Figure 3-9b. Optimum resonator-transducer separation 
distance (Ino) for maximum power gain versus frequency. 
Po = 1 atm. 
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E. ANALYSES WITH TRANSDUCER ka = CONSTANT AND Pq = 50 
ATM 

Two major assumptions that have been made throughout the 
analysis up this point are now modified. The first is that the ambient 
pressure is one atmosphere and the second is that the radius of the 
transducer equals that of the resonators (kao = kan). In the following 
analysis, an ambient pressure of 50 atmospheres (corresponding to a 
depth of 500 meters) will be taken as a more realistic value for an 
operational system. Since the internal acoustic impedance of a bubble 
is a function of depth, the value for kar of a single bubble at resonance 
will shift from 0.01379 (Pq = 1 atm) to 0.09754 (Pq = 50 atm). Also, 
the value of kao will no longer be considered equal to that of the 
resonators, but will be held constant at various specific values. 

In the previous sections it was shown that little is gained when the 
number of sympathetic resonators is greater than about four; hence 
the following analysis will be restricted to the case N = 4. Conducting 
similar calculations as before, contour plots and surface plots have 
been generated for the cases N = 4, kao = 0.05, 0.10 and 0,20, for 
both the conical configuration (Figures 3- 10a through 3-lOf) and the 
planar configuration (Figures 3- 11a through 3-1 If). As in the previous 
analyses where the ambient pressure was taken to be one atmosphere, 
all cases display a local maximum of Gr of slightly less than 2 for ka s 
kar = 0.09754 and kl = n. These mciximum values are tabulated in 
Table 3-3. Also tabulated, though not plotted, are the cases kao = 
0.025 and 0.150. The first entry in Table 3-3 for either configuration 
was calculated by holding kao equal to kan, as in previous calculations, 
vice holding kao constant. A review of Table 3-3 and Figures 3-11 and 
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12 demonstrates that the value of kao has little effect on the achievable 
power gain. 



configuration 


kao 


Gr 


kamax 


^max 


magOJSri) 


arg(Ui,'‘) 


planar 


kUn 


1.87 


0.09987 


2.82 


0.649 


90.0° 




.025 


1.87 


0.09987 


2.82 


0.652 


90.0° 




.050 


1.87 


0.09987 


2.82 


0.652 


90.0° 




.100 


1.87 


0.09987 


2.82 


0.649 


90.0° 




.150 


1.87 


0.09987 


2.82 


0.645 


90.0° 




.200 


1.87 


0.09998 


2.82 


0.639 


90.0° 


conical 


ka„ 


1.78 


0.09990 


3.08 


0.601 


90.0° 




.025 


1.78 


0.09990 


3.08 


0.603 


90.0° 




.050 


1.78 


0.09990 


3.08 


0.603 


90.0° 




.100 


1.78 


0.09990 


3.08 


0.601 


90.0° 




.150 


1.78 


0.09990 


3.08 


0.597 


90.0° 




.200 


1.78 


0.09990 


3.08 


0.592 


90.0° 



Table 3-3. Values of Gr for various values of kao ranging 
from 0.025 to 0.2 and the parameters associated with the 
local maximum at kl = n, ka = kar = 0.09754. Conical and 
planar configurations. Pq = 50 atm. 
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Figure 3-10a,b. Surface and contour plots of Gr as a 
function of ka and klno for N = 4. Planar configuration. Pq 
= 50 atm. kao = 0.05. 
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c. 




Figure 3-10c,d. Surface and contour plots of Gr as a 
function of ka and klnO for N = 4. Planar configuration. Pq 
= 50 atm. kao = 0.10. 



48 



e. 




Figure 3-10e,f. Surface and contour plots of Gr as a 
function of ka and klno for N = 4. Planar configuration. Pq 
= 50 atm. kao = 0.20. 
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a. 




Figure 3-lla,b. Surface and contour plots of Gr as a 
function of ka and kl^o for N = 4. Conical configuration. 
Pq = 50 atm. kao = 0.05. 
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Figure 3-llc,d. Surface and contour plots of Gr as a 
function of ka and klnO for N = 4. Conical configuration. 
Pq = 50 atm. kao = 0.10. 
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Figure 3-lle,f. Surface and contour plots of Gr as a 
function of ka and klno for N = 4. Conical configuration. 
Pq = 50 atm. kao = 0.20. 
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E.l. BEAM PATTERNS, Pq = 50 ATM 

Since there is no signifieant variation in gain with the value of kaQ, 
beam patterns were only developed for the case kag = kaj-, and are 
shown in Figures 3-12 through 3-14. These plots show the maximum 
SPL in the positive z direetion for the conical configuration to be 1.5 
dB greater than for the planar case and 10.5 dB greater than for the 
transducer in the absence of the resonators. 



z axis 




Figure 3-12. Beam patterns In the y,z plane for tlie planar 
and conical configurations. Pg = 50 atm. Scale in decibels 
(dB). 
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z axis 



o 




Figure 3-13. Beam patterns in the y,z and z,0 = 45° 
planes. Planar configuration. Pq = 50 atm. Scale in 
decibels (dB). 
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o 




Figure 3-14. Beam patterns in the y.z and z,0 = 45° 
planes. Conical configuration. Pq = 50 atm. Scale in 
decibels (dB). 
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E.2. OPTIMUM CONFIGURATION, Pq = 50 ATM 

The actual dimensions of the resonator radius a and the separation 
distances 1^0. for a system operated at an ambient pressure of 50 
atmospheres, were calculated using Equation 3-5 using the optimum 
values ka = kar = 0.09754 and klno s n. Since the values of Ino are not 
depth-dependent. Figure 3-9b still applies. Figure 3-15 displays the 
relation between the optimum resonator radius and frequency for an 
ambient pressure of 50 atmospheres. 
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Figure 3-15. Optimum resonator diameter (a) at 
resonance verses frequency. Pq = 50 atm. 
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IV. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 



A. SUMMARY 

The problem of multiple-scattering of an incident plane-wave from 
a system of compact resonant-scatterers was solved using both 
Tolstoy’s method and network analysis. It was shown that these two 
methods of analysis are equivalent and that both predict the existence 
of "quasiresonance". Network analysis was then applied to a circular 
array of scatterers, where the incident plane-wave was replaced by 
spherical waves radiated from an active transducer, located on the 
center axis. The acoustic advantage of this system was quantified by 
the gain in the radiation resistance seen by the transducer in the 
presence of these scatterers compared to a lone transducer. Two 
configurations were analyzed: (1) the transducer in the plane of the 

resonators, and (2) the transducer displaced by one-quarter 
wavelength from the plane. Both configurations showed a gain in 
radiation resistance of almost a factor of two and an improved 
directivity over that of a single transducer. 

B. CONCLUSIONS 

The numerical analysis of a transducer surrounded by a circular 
array of sympathetic resonators indicates that a gain in radiation 
resistance of about two is obtainable compared to a lone transducer. 
The physical dimensions of such a system operating at an ambient 
pressure of 50 atmospheres are displayed in Table 4-1 for several 
frequencies of interest. 
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Frequency 

(Hz) 


Resonator Radius (a) 
(m) 


System Radius (R 
(m) 

conical planar 


10 


0.329 


64.9 


74.9 


100 


0.0329 


6.49 


7.49 


1000 


0.00329 


0.649 


0.749 



Table 4-1. Dimensions of planar and conical 
configurations for various frequencies. Pq = 50 atm. R is 
the radius of the array. 



Although the calculations presented here were for spherical 
bubbles, the results also apply to more realistic devices (e.g. thin 
shelled cavities, balloons, etc) as long as they are compact (ka « 1). 
M. Strasberg [Ref. 6] showed that the shape of a resonant body has 
little effect on its fundamental resonance frequency. For example, he 
calculated that an oblate spheroid with a ratio of lengths of major to 
minor axes of four, the fundamental resonance frequency differs by 
only a factor of 1.08 from that of a sphere with the same volume. 
Therefore, any deformation of the resonators from a true spherical 
shape, for example due to forces such as buoyancy, should have little 
affect on system performance. 



C. RECOMMENDATIONS 

The results of the research reported in this thesis are very 
encouraging. It is recommended that follow-on work include the 
design, construction, and testing of an array to verify the theoretical 
predictions presented. 

Theoretical analysis remains to be conducted in at least four areas. 
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First, in all of the calculations performed, the sympathetic resonators 
were assumed to be gas-filled air bubbles executing radial oscillations 
only. A more realistic model of a sympathetic resonator may be 
required, based upon the properties of a practical device. It is 
recommended this be investigated. Second, Tolstoy has shown that 
the scattering of an incident plane-wave from a linear array of resonant 
scatterers exhibits strong directionality. Perhaps a linear array of 
sympathetic resonators offers better performance gain than the 
circular sirrays considered in this research. It is recommended that 
such a configuration be investigated. Third, it was shown that in the 
optimum conical configuration the ring of sympathetic resonators lies 
on a cone with a apex angle of approximately 120 degrees, 
independent of frequency. This suggests that a broad-band, constant 
directivity system could be constructed using a series of concentric 
rings of resonators, each tuned to a slightly different frequency and 
placed in its appropriate location along the surface of such a cone, as 
indicated in Figure 4-1. It is recommended that the properties of 
such a system be investigated. Fourth, in the calculations presented, 
the sympathetic resonators were all assumed to be identical, and their 
locations were assumed to be precisely known, as well as that of the 
transducer. It is recommended that the effects of resonator 
nonuniformity and imprecise location be investigated. 
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Transducer 




Figure 4-1. Multiple-ring conical configuration allowing an 
expanded frequency response. Two concentric rings of 
resonators of radii a and a’, separated from the transducer 
by distances of 1 and 1’, respectively. 
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LIST OF SYMBOLS 



Radius of element n. 

c Speed of sound in water. 

Gr The ratio of the acoustic radiation resistance of a transducer 
with sympathetic resonators present to the acoustic radiation 
resistance of an identical transducer with no resonators 
present. 

k Wave number, co/c. 

kaj- The resonance value of ka for a single air- filled bubble. 

^nm Separation distance between element n and m. 

pfi Incident free field acoustic pressure at element n due to 

external field only. 

Pn Surface pressure at element n. 

p&c Scattered pressure at the surface of resonator n, in the 
presence of other resonators. 

Pifngle Scattered pressure at the surface of a single resonator, i.e. in 
the absence of other resonators. 



r 

z 



A 

B 

^nm 



Distance from scatterer center to a field point. 

Specific acoustic impedance, pressure divided by surface 
normal velocity. 



Single-scattering amplification coefficient of element n. 



Multiple-scattering amplification coefficient of element n. 



Ratio of acoustic impedances, defined as; Cnm 



-7a r 
^nm 
-yar 
^mm 



Ratio of acoustic impedances, defined as: Dn 



yar . 
^nn ^ 



Zg 



'7ar 
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Fa Scattered acoustic pressure from a single scatterer, 

normalized to the incident acoustic field. 

Scattered acoustic pressure from a scatterer, interacting with 
other scatterers, normedized to the incident acoustic field 

I Electrical current. 

Pq Ambient pressure of environment. 

S„ Surface area of element n. 

T Transduction coefficient. 

Un Volume velocity of element n. 

Volume velocity of element n, normalized to transducer 
volume velocity. 

Z Acoustic impedance, pressure divided by volume velocity. 

Zn Acoustic impedance of element n. 

Znn Acoustic radiation impedance seen by element n with all 
other elements blocked. 

Znm Acoustic transfer impedance from element m to element n. 

Zfot The total acoustic radiation impedance seen by transducer in 
the presence of sympathetic resonators. 

Z^o Transducer open circuit (1=0) acoustical impedance. 

^nm Kronecker delta function, where: 6nm = | ^ . 

I 0; m ^ n I 

p Density of water. 

p Ratio of the magnitudes of the scattered acoustic pressure at 

the surface of a resonant scatterer to the incident acoustic 
pressure in the presence of other scatterers. 
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|io Ratio of the magnitudes of the scattered acoustic pressure at 

the surface of a single resonant scattererto the incident 
acoustic pressure. 

Y Ratio of specific heats (air). 

(0 Angular frequency of plane wave. 

(Oq Angular frequency of scatterer resonemce. 

^ Defined as: ^ ^ . 

^ ka 

IXjjl A Nxl column matrix, n = 1...N. 

[Xnml A NxN square matrix, n,m = 1...N. 
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